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Abstract 

We investigate four-dimensional compact U(l) lattice gauge theory with 
a monopole term added to the Wilson action. First we consider the phase 
structure at negative /?, revealing some properties of a third phase region there, 
in particular the existence of a number of different states. Then our present 
studies concentrate on larger values of the monopole coupling A where the 
confinement-Coulomb phase transition turns out to become of second order. 
Performing a finite-size analysis we find that the critical exponent is close 
to, however, different from the gaussian value and that in the range considered 
v increases somewhat with A. 



I. INTRODUCTION 



The phase transition in 4-dimensional compact U(l) lattice gauge theory with the 
Wilson action because of the occurrence of a gap in the energy histogram is believed 
to be of first order. Actually this is to be studied in more detail by a finite-size 
analysis. Recently critical exponents have been determined in a higher-statistics 
study There the critical exponent v has been found to decrease towards \ with 
increasing lattice size, i.e. towards the value characteristic of a first-order transition. 

If the Wilson action is extended by a double charge term with coupling 7 the first 
order transition weakens with decreasing 7 and has been conjectured to become of 
second order at some negative 7 [0] , which, however, so far has not been confirmed. 
In Ref. instead of the usual periodic boundary conditions a spherelike lattice 
(i.e. the surface of a 5-cube which is homeomorphic to a 4-sphere) has been used. 
The disappearance of the gap at 7 = 0, -0.2, -0.5 observed there has been interpreted 
by the authors as an earlier start of the second order region and the relatively low 
value of V there as related to a nongaussian fixed point. 

In a more recent investigation of the action with double charge term the gap 
has been shown for 7 down to —0.4 to reappear on larger spherelike lattices. The 
absence of a gap on smaller lattices has been attributed to larger finite-size effects 
in the spherical geometry. These observations are just what was to be expected 
according to studies of the influence of different boundary conditions where it 
turned out that inhomogeneities weaken the transition and that the gap reappears 
for sufficiently large lattice sizes. In addition in Ref. 0] for 7 between +0.2 and 
—0.4 the critical exponent v has been found to decrease towards | with increasing 
lattice size for toroidal as well as for spherical geometry. Further, in some cases 
stabilization of the latent heat has been observed. 

Thus, for the above actions there are now rather strong indications that, at least 
in the region which has been investigated so far, the transition is of first order. Of 
course, it remains desirable to check this on still larger lattices. While this is hardly 
possible with conventional methods, it can be done using the dynamical-parameter 
algorithm developed in @, 0| . 

After the mentioned evidence of first order the question arises where one can 
find a second-order phase transition in the U(l) gauge system. This leads to the 
modification of the Wilson action where a monopole term is added , which reads 

S = ^Y. (l-cose^,,,) + A^|M^,,| (1.1) 
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with Mp^^ = epaf^u{Qtiu,x+a - Q,iu,x)/4:n where the physical flux Q^^^^ G [-7r,7r) is 
related to the plaquette angle 6^,,,^ G {-Att,Att) by Q^^^^ = Q^^^^ + inrip^^^ 0. 

Using the action ( |1 . 1| ) from the energy distribution it is seen that the gap gets 
smaller with increasing monople coupling A 0], which indicates that the flrst- 
order transition gets weaker and for sufficiently large A becomes of second order. 
The latter has been corroborated by a flnite-size analysis [jlO| which has shown that 
at A = 0.9 the critical exponent is already characteristic of a second-order transition. 

In addition to showing a second-order phase transition for sufficiently large A, the 



action (IJ.) is particularly attractive in view of the fact [0, 11] that the confinement 
phase and the Coulomb phase are unambiguously characterized by the presence or 
absence, respectively, of an infinite network of monopole current lines, where "infi- 
nite" on finite lattices is to be defined in accordance with the boundary conditions 
0. Since the probability Pnct to find an infinite network takes the values 1 and 0, 
respectively, it is very efficient to discriminate between those phases. For the finite 
lattices with periodic boundary conditions we are considering here, "infinite" means 
"topologically nontrivial in all directions". While for loops this characterization is 
straightforward, to determine the topological properties of the occurring networks 
it is necessary to perform a more elaborate analysis based on homotopy preserving 
mappings [H, |lT|. 

In |jlO| the location of the transition from the confinement phase to the Coulomb 
phase has been determined as a function of A. It has been found that the related 
transition line continues to negative (5. On the other hand, in Refs. |T2|, |T3| a further 
transition has been reported at /3 = — 1 for A = and at /? = — 0.7 for A = oo. This 
suggests to check in more detail what happens at negative (3. 

In the second-order region mentioned above an important question is wether one 
has universal critical properties there. The energy distribution indicates that this 
region starts at some finite A above 0.7 [§, |^. For A = 0.9 the critical exponent v is 
known to be characteristic of a second-order transition [1^. It appears important 
now to get more information on this region. In particular finite-size analyses at 
different values of A and for a variety of observables are desirable. 

In the present work we have performed Monte Carlo simulations to address the 
indicated questions. Simulation runs at various values of A and {3 have been done 
to get an overview of the situation at negative (3. The emerging picture of this is 
discussed in Sect. II. The emphasis of our investigations has then been on higher- 
statistics simulations at A = 1.1 and A = 0.8 in the critical region of the confinement- 
Coulomb phase transition, which have been evaluated by finite-size analyses. The 
results of this are presented in Sect. III. 
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II. PHASE REGIONS 



In the location (3cr of the transition from the confinement phase to the 
Coulomb phase has been determined by the maximum of the specific heat Cmax 
up to A = 1.3. Because the peak of the specific heat has been found to decrease 
strongly with A, at larger values of A only Pnet has been used to locate the transition. 
In these investigations it has turned out that above A = 1.2 the associated values of 
Per become negative. 

The Wilson action has the symmetry (3 — >• — /3, — > — Q. At A = it gives 



rise to a phase transition at P = —1 in addition to the one at /9 = 1 |12]. For A 7^ 



this symmetry is violated by the monopole term. At A = cxo only the transition at 



negative /3 persists and occurs at about (3 = —0.7 |T^, [T3|. 

Here we have checked the occurrence of such transition at negative (3 also at 
intermediate values of A by determining Cmax- Our observations indicate a transition 
line extending from (A, (3) = (0, —1) to (00, —0.7). In the energy distribution on the 
8^ lattice we have seen a double-peak structure at A = 0.5 and indications of such 
structures (though not resolved with present accuracy) at other A values. Together 



with the observations at A = and at A = 00 in Refs. [|^, [T3| this hints at first 
order along that transition line. Of course, it remains to confirm this with higher 
statistics and on larger lattices. 

We have also investigated some properties of the system below the new transition 
line which will be discussed in the following. Since so far we have no indication of a 
further subdivision of this region, but rather find similar properties throughout it, 
we consider it here (at least as a working hypothesis) as a third phase. 

Figure 1 gives an overview of the phase regions as they are according to our present 



knowledge. The line separating confinement and Coulomb phases obtained in |]10 
with Cmax has been supplemented by a point at A = 1.35 determined here and by 
the point (1.4, 0.52) found there only using Pnet- It is seen that this transition line 
hits the boundary of the third phase at approximately (1.43, —0.71). 

While Pnet has provided an unambiguous criterion in the confinement phase and 
in the Coulomb phase, this is no longer so in the third phase. At fixed (A,/3) in 
this region values as well as 1 occur for Pnet- The observations at (10, —1000) in 
HlOl appear to be related to this. That properties of monopole structures are no 



^The transformation Uo —Uo of the plaquette fields can be realized by transforming the link 
angles as Q^^x ^^fix^fix where the sign choice is according to Qfj^x being in the intervals [— tt, 0) 
and [0, tt), respectively, and where A^^; = — 1 if xi + 2:2 + • • • + s^^^-i is odd and A^^^ = +1 if it is 
even or if /i = 1. 
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longer characteristic at sufficiently negative /3 can be understood by noting that the 
respective quantities are not invariant under the transformation (5 — > — /3, [/□ — > 
-Uu. 

As a characteristic feature of the third phase we have found that a number of 
different states exists there between which transitions in the simulations are strongly 
suppressed. We have observed this phenomenon at various negative (3 in the A 
range from to 2.5 . Typical examples of time histories of the average plaquette 
e = (1/6L*^) X]^>!/,x(l ~ cosO^jy^a;) are given in Figure 2, with the results of seven 
simulation runs, exhibiting different states and some transitions from lower to higher 
ones. It turns out that there is no correlation between the state reached and the 
type of start (hot or cold) of a simulation run. 

Comparing time histories of e at different (A, /5) the splittings of the states appear 
very similar. Closer inspection shows that at least nine states occur. The number of 
transitions between states increases with A. They have so far always been observed 
occurring from lower to higher plaquette values. The widths in the time histories 
show httle dependence on A while they strongly decreases if (3 gets more negative 
(width x/5 being roughly constant). Thus at smaller negative j3 it gets increasingly 
difficult to resolve splittings. 

The average e found show not much dependence on A, however, they increase 
as (3 gets more negative. The e values, being somewhat below 2, indicate that 
at sufficiently negative 13 the average of cosQ^u,x gets close to —1. That negative 
cos Q^y^x occur at negative [3 and positive cos Q^v^x at positive [3 in view of the 
symmetry j3 —f3, Uu — ^ —Uu of the Wilson action is conceivable. 

Though at negative j3 the monopole density p = (1/4L'^) Y,p,x \ ^pA appears of less 
interest, it should be mentioned that the states arc also seen in its time histories. 
The respective time histories are correlated with large values of p corresponding to 
low values of e. The widths for p decrease with A. They are rather large so that, 
except at the largest values of A, the resolution of the states is very poor. The value 
of p decreases with A. 

The origin of the observed different states is not yet clear. Similar observations 
are made in simulations of spin glasses and of frustrated systems, and also in finite- 
temperature SU(N) gauge theory with states related to spontaneous Z{N) breaking. 
In any case the third phase appears to have more complicated properties which 
deserve further (though computationally demanding) studies elsewhere. 
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III. CRITICAL PROPERTIES 

Monte Carlo simulations with about 10^ sweeps have been performed at each of up 
to 20 values of (3 in the critical region of the confinement- Coulomb phase transition 
at A = 1.1 and at A = 0.8 and for each lattize size considered. Multihistogram 
techniques |1^ have been applied in the evaluation of the data. The errors have 



been estimated by Jackknife methods. In the finite-size analyses in addition to the 
specific heat and the Challa-Landau-Binder (CLB) cumulant [|r3] complex zeros of 
the partition function [|T6|, in particular the Fisher zero zq closest to the [3 axis p^ , 
have been used. 



The specific heat is 

where E = J2fi>u,x{^ ~ cos Q^u,x)- For d = 4 its maximum is expected to behave as 



C = -^,{{E - {E)r) (3.1) 



Cmax ~ (3.2) 

if the phase transition is of first order and as 

Cmax ~ L~ (3.3) 

if it is of second order, where a is the critical exponent of the specific heat and u 
the critical exponent of the correlation length. 

In Figure 3 we present the results for Cmax obtained on lattices with L = 6, 8, 
10, 12, 16 for A = 0.8, 0.9, 1.1. They include data from simulations of the present 
investigation (A = 1.1 and A = 0.8) and ones from simulations of Ref. [|10| (A = 0.9). 
From Figure 3 it is already obvious that ^ as it occurs in ( p.3|) decreases with A, 
which by the hyperscaling relation a = 2 — d u means that u increases with A. 

The fits to the data presented in Figure 3 give the values for ^ shown in Table I. 
They are clearly very far from 4 and thus the transition cannot be expected to be 
of first order. Using the hyperscaling relation a = 2 — du the values for u listed 
in Table II are obtained. They are different from the value | of the gaussian case, 
however, rather close to it. Thus in any case to conclude on second order appears 
quite safe. 

Similar results are obtained for the minimum of the CLB cumulant 



3^ {E^y 
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and for the imaginary part of the closest Fisher zero zq. For these quantities finite- 
size scahng predicts the behaviors 



Im(zo) ~ " , (3.5) 

f^CLB.min ~ L~~^ . (3.6) 

The resuhs of the respective fits are also listed in Tables I and II. It is seen that 
the values obtained from different quantities roughly agree, with some systematic 
deviations beyond the given statistical errors. In any case it is obvious that u 
increases somewhat with A and that it is not far from |. 

The related critical (3cr{L) (i.e. the extrema positions of C and f/cLB and the real 
part of zq) are given in Table III . They illustrate the dependence on A and lattice 
size. The critical (3cr{L) are expected to behave as 

(3,,{L) = (3^^ + aL-"- . (3.7) 

Using the values of u in Table II and the data for (3cr{L) in Table III, the numbers 
Pct{oo) and a in ( p^) have been calculated for A =1.1. In Table IV they are compared 



with those obtained for A=0.9 in |T0[. The errors given are again only statistical 
ones. 

The observed increase of u with A could indicate a nonuniversal behavior with a 
maximal value (which could be even |) reached at the boundary to the third phase. 
Another possiblity is that, because the observations are on finite lattices and the 
range of A considered is not too far from the first-order region, the starting point of 
the second-order region is not yet sharp. Then this effect should disappear on much 
larger lattices. 

To get hints on the dependence on lattice size effective u (i.e. a sequence of ones 
related to two neighboring lattice sizes) have also been considered. Though there 
was a slight tendency in support of the above view of a finite-size effect, it turns out 
that larger lattices and higher statistics (which unfortunately are computationally 
very expensive) would be needed for firm conclusions. 

If for very large lattices u gets independent of A for the A above a certain finite 
value and below the start of the third phase, the question arises what the ultimate 
u would be. Because the numbers for u in Table II come rather close to | this could 
be simply the gaussian value. If first order is confirmed along the boundary to the 
third phase, an interesting apect is that this transition line with z/ = | is met by 
the second-order line with u close to or equal to |. 
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TABLE I 



Critical exponents ^ from Cy, Uclb- 



A 


Cv 


^CLB 


0.8 


0.616(22) 


0.756(23) 


0.9 


0.485(35) 




1.1 


0.284(12) 


0.391(12) 



TABLE II 



Critical exponents u from Im(2;o), Cy, C^clb- 



A 


Im(2;o) 


Cv 


t^CLB 


0.8 


0.404(5) 


0.433(2) 


0.421(3) 


0.9 




0.446(5) 




1.1 


0.421(8) 


0.467(2) 


0.455(2) 
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TABLE III 



Critical (3cr from zq, Cy, ^/clb- 



A 


L 


Re(2;o) 


Cv 


%CL 


0.8 


6 


0.4637 (3) 


0.4630 (3) 


0.4647 (2) 




8 


0.4783 (2) 


0.4779 (2) 


0.4783 (4) 




10 


0.4844 (2) 


0.4841 (2) 


0.4843 (2) 


1.1 


6 


0.1409 (10) 


0.1387(5) 


0.1431 (6) 




8 


0.1636 (10) 


0.1618(5) 


0.1630 (6) 




10 


0.1742 (4) 


0.1731(5) 


0.1736 (6) 




12 


0.1775 (8) 


0.1782(8) 


0.1784 (8) 




16 


0.1839 (1) 


0.1837(1) 


0.1838 (1) 



TABLE IV 



Fit data (3^ and a for zq^ Cv-, ^/clb- 



A 






Cv 


^CLB 


0.9 


A^cr 




0.4059(5) 






a 




-1.99(6) 




1.1 


f^cr 


0.1888 (6) 


0.1902 (5) 


0.1896 (5) 




a 


-3.14(13) 


-2.41 (7) 


-2.49 (8) 
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Figure captions 



Fig. 1. Location of phase transition points Per on 8^ lattice as function of A 

between confinement and Coulomb phases (circles from Cmaxj square from 
and to third phase (diamonds). Curves drawn to guide the eye. 

Fig. 2. Typical time histories of the average plaquette e at negative P 
obtained for seven different simulation runs, 
shown for (A,/3) = (2, -20) and 8^ lattice. 

Fig. 3. Maximum of specific heat Cmax as function of lattice size L 
for A = 0.8, 0.9, and 1.1 at transition points Per between 
confinement and Coulomb phases. 
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Figure 2: 
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Figure 3: 
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